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Abstract 

Steady-state  models  are  developed  to  describe  an  experimental  methane  fuel  processor  that  is  intended  to  provide  hydrogen  for  a  fuel 
cell  system  for  power  generation  (2-3  kW).  First-principle  reactor  models  are  constructed  to  describe  a  series  of  reactions,  i.e.,  steam 
and  autothermal  reforming  (SR/ATR),  high-  and  low-temperature  water-gas  shift  (HTS/LTS)  reactions  and  preferential  oxidation  (PROX) 
reactions,  at  different  sectors  of  the  reactor  system  for  methane  reforming  as  well  as  gas  cleaning.  The  pre-exponential  factors  of  the  rate 
constants  are  adjusted  to  fit  the  experimental  data  and  the  resultant  reactor  model  provides  a  reasonably  good  description  of  steady- state 
behaviour.  Next,  sensitivity  analyses  are  performed  to  locate  the  optimum  operating  point  of  the  fuel  processor.  The  objective  function  of  the 
optimization  is  fuel  processor  efficiency.  The  dominating  optimization  variables  include:  the  ratios  of  water  and  oxygen  to  the  hydrocarbon 
feed  to  the  autothermal  reforming  reactor  and  the  inlet  temperature  of  the  reactor.  The  results  indicate  that  further  improvement  in  fuel 
processor  efficiency  can  be  made  with  a  reliable  process  model. 

©  2005  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Polymer  electrolyte  fuel  cell;  Fuel  processor;  Optimization;  Anthermal  reforming;  Steam  reforming;  Preferential  oxidation 


1.  Introduction 

Fuel  cell  systems  have  received  a  great  deal  of  atten¬ 
tion  in  recent  years  because  of  their  competitive  efficiency 
and  reduced  emissions  [1].  Fuel  cell  systems  can  be  used  in 
stationary  or  mobile  applications.  The  polymer  electrolyte 
membrane  fuel  cell  (PEMFC)  is  attracting  much  attention; 
it  uses  either  pure  hydrogen  or  a  hydrogen-rich  gas  mixture 
for  power  generation.  Because  of  the  difficulty  associated 
with  hydrogen  storage  [2],  a  reformer  is  often  employed  to 
convert  hydrocarbons  to  hydrogen-rich  synthesis  gas  (‘syn¬ 
gas’).  The  syngas  leaving  out  of  the  reformer  usually  con¬ 
tains  hydrogen,  carbon  dioxide,  water  vapour,  unconverted 
fuel  and  carbon  monoxide  (CO).  Because  of  stringent  limi¬ 
tation  on  the  carbon  monoxide  concentration  in  the  PEMFC 
[2],  a  series  of  CO-removing  steps,  such  as  the  water-gas 
shift  reaction  and  the  preferential  oxidation  reaction,  are 
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adopted  to  keep  CO  concentration  below  100  ppm  before  en¬ 
tering  the  fuel  cell.  A  reformer  and  several  CO-removing 
reactors  constitute  a  fuel  processor.  The  fuel  for  the  fuel 
processor  will  vary  with  different  applications.  It  may  be  a 
liquid  fuel  such  as  methanol,  gasoline  and  diesel  for  trans¬ 
portation  systems,  or  natural  gas  or  propane  for  stationary 
systems  [1].  Options  for  converting  fuels  to  hydrogen-rich 
syngas  in  the  reformer  include:  steam  reforming  (SR)  which 
uses  fuel  and  water  as  feed;  partial  oxidation  (POX)  which 
uses  fuel  and  oxygen  as  feed;  autothermal  reforming  (ATR) 
which  uses  fuel,  water  and  oxygen  as  feed,  i.e.,  a  combina¬ 
tion  of  the  endothermic  SR  reaction  and  the  exothermic  POX 
reaction  [3]. 

Reviews  on  various  aspects  of  fuel  processing  can  be  found 
in  [4-7] .  Research  emphasis  is  still  placed  upon  the  reforming 
of  hydrocarbon  fuels.  Reaction  kinetics  and  reactor  models 
have  also  been  studied  in  [8-10].  A  comparison  of  the  ef¬ 
ficiency  for  different  fuels  has  been  explored  in  [1],  while 
simulation  and  reactor  performance  analyses  have  been  ex¬ 
amined  in  [3,10-12].  Relationships  between  thermodynamic 
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properties  and  the  expected  performance  of  fuel  processors 
have  been  presented  in  [13,14]. 

This  study  presents  a  systematic  approach  towards  the 
development  and  validation  of  a  model  for  an  experimen¬ 
tal  fuel  process.  First,  the  necessary  experimental  steps  are 
performed.  Then,  a  steady-state  model  is  established  and  the 
kinetics  parameters  are  adjusted  to  provide  a  description  of 
the  system  behaviour.  Finally,  a  sensitivity  analysis  is  per¬ 
formed  to  determine  the  optimum  operating  conditions  for 
improved  efficiency. 

2.  Experimental 

An  experimental  fuel  processor  was  designed  and  installed 
in  the  facility  of  the  Union  Chemical  Laboratory  (UCL)  of 
the  Industrial  Technology  Research  Institute  (ITRI).  The  en¬ 


tire  fuel  processor  set-up  is  composed  of  several  reactors, 
heat-exchangers  and  a  burner.  A  detailed  drawing  of  the  sys¬ 
tem  with  methane,  water  and  oxygen  as  feed  is  presented 
in  Fig.  1.  The  various  feeds  enter  the  system  at  room  tem¬ 
perature  and  their  individual  flow  rates  can  be  adjusted.  The 
mixed  feed  first  passes  through  a  heat-exchanger,  the  tem¬ 
perature  of  which  increases  slightly.  Next,  the  feed  passes 
through  a  pipe  on  the  top  of  the  reformer,  which  is  coiled 
down  in  the  outside,  and  then  enters  the  reformer.  A  burner 
is  placed  under  the  reformer.  It  also  uses  methane  as  fuel 
and  supplies  heat  directly  to  the  reformer  for  pre-heating  and 
reaction  purposes.  The  reformer  is  made  of  stainless  steel 
304  that  has  a  good  heat-transfer  capability.  A  honeycomb 
catalyst  carrier  is  placed  inside  the  reformer  and  the  catalyst 
is  coated  on  a  ceramic  carrier.  The  catalyst  was  developed 
by  UCL  of  ITRI  and  is  a  Ru/Ce02-Zr02  catalyst.  The  re¬ 
former  and  burner  are  wrapped  together  in  insulation  mate- 


Fig.  1.  Detail  drawing  of  experimental  fuel  processor:  Tl,  temperature  of  burner;  T2,  reformer  inlet  temperature;  T3,  reformer  outlet  temperature;  T4,  HTS1 
inlet  temperature;  T5,  HTS2  inlet  temperature;  T6,  LTS  inlet  temperaturel;  T7,  PROX1  inlet  temperature;  T8,  PROX2  inlet  temperature. 
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rials  to  prevent  heat  loss.  The  reactor  effluent  exits  out  from 
the  top  of  the  reformer  and  is  further  processed  in  a  series 
of  water-gas  shift  reactors  and  preferential  oxidation  reac¬ 
tors.  The  hot  reformer  effluent  is  cooled  down  by  means  of  a 
heat-exchanger  and  direct  water  injection  before  entering  the 
water-gas  shift  reactors.  The  temperatures  of  the  water-gas 
shift  reactors  are  arranged  in  decreasing  order,  i.e.,  water  in¬ 
jection  between  reactors  is  devised  to  cool  down  the  syn¬ 
gas.  There  are  two  high- temperature  water-gas  shift  reactors 
(HTS1  and  HTS2)  and  one  low-temperature  water-gas  shift 
reactor  (LTS).  The  operating  temperature  is  between  300  and 
400  °C  for  the  HTSs,  which  contain  a  Pt/mixed- oxide  cata¬ 
lyst  that  has  been  developed  by  UCL  of  ITRI  and  is  coated  on 
honeycomb  carriers.  The  operating  temperature  of  the  LTS  is 
between  190  and  250  °C  and  a  commercial  Cu-Zn0/Al203 
catalyst  is  used.  Oxygen  is  introduced  directly  before  the  syn¬ 
gas  enters  the  preferential  oxidation  reactors  (PROX).  The 
reactions  are  promoted  by  a  catalyst  that  has  been  prepared 
by  UCL  of  ITRI  and  contains  3  wt.%  Pt. 

The  objectives  of  the  experimental  set-up  are: 

(i)  to  test  the  reactivity  and  property  of  the  catalysts  devel¬ 
oped  by  UCL  of  ITRI; 

(ii)  to  obtain  operating  parameters  to  achieve  the  desired 
hydrogen  flow  rate  and  carbon  monoxide  concentration; 

(iii)  to  compare  the  efficiencies  of  the  different  reaction  path¬ 
ways  (auto thermal  reforming  and  steam  reforming). 

The  experiment  is  operated  in  three  stages:  (i)  heat-up, 
(ii)  autothermal  reforming,  (iii)  steam  reforming.  In  the  first 
5  min,  methane  and  air  (oxygen)  is  fed  into  the  reformer  to 
carry  out  the  combustion  reaction  and  a  large  quantity  of 
fuel  is  added  to  the  burner  to  heat  up  the  catalyst  and  the 
monolith  support  of  the  reformer.  The  purpose  of  this  stage 
is  to  heat  the  reformer  to  the  reaction  temperature.  The  next 
step  is  to  add  water  and  the  reformer  follows  the  autothermal 
reforming  (ATR)  pathway.  A  suitable  amount  of  fuel  is  added 
to  maintain  the  reformer  inlet  temperature  at  about  667  °C. 
Finally,  the  air  feed  is  stopped  after  3.5  h  with  only  methane 
and  water  as  feeds.  This  implies  that  because  of  a  higher 
operating  temperature  as  well  as  strong  exothermic  reaction, 
the  amount  of  fuel  fed  to  the  burner  is  larger  than  that  required 
in  the  pathway.  The  results  of  typical  experimental  runs  are 
given  in  Figs.  2-4.  The  fuel  flow  rates  and  temperatures  at  the 
inlet  and  the  outlet  of  the  reformer  at  the  ATR  and  SR  stages 
are  shown  in  Fig.  2 A  and  B,  respectively.  The  temperature 
profiles  for  different  sectors  of  the  fuel  processor  are  given  in 
Fig.  3,  and  the  concentration  profiles  for  CO,  hydrogen  and 
methane  conversion  in  Fig.  4. 


3.  Steady-state  modelling  of  fuel  processor 

3.1.  Description  of  system 
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Fig.  2.  Feed  CH4  flow  rate  (reformer  and  burner),  reformer  inlet  and  outlet 
temperature  at:  (A)  start-up  and  ATR  pathway  stages  and  (B)  SR  pathway 
stage. 

and  includes  a  reformer,  a  heat-exchanger  and  a  burner  in  a 
heat-integrated  framework,  as  shown  in  Fig.  6.  The  second 
part  is  the  gas-cleaning  unit,  which  is  composed  of  reactors 
connected  in  series.  Three  of  them  carry  out  the  water-gas 
shift  reaction  and  two  perform  the  preferential  oxidation  re¬ 
action.  These  reactors  are  lumped  together  in  the  modelling 
phase. 
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The  general  structure  of  the  experimental  fuel  processor 
consists  of  two  major  parts.  The  first  part  is  the  reforming  unit 


Fig.  3.  Inlet  temperatures  of  HTS1,  HTS2,  LTS,  PROX1  and  PROX2  at 
start-up  and  ATR  pathway  stage. 
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Fig.  4.  CO  concentration,  methane  conversion  and  H2  concentration  at  re¬ 
former  outlet  and  various  reactor  outlets  in  gas-cleaning  unit. 


3.2.  Reactor  modelling  and  corresponding  reaction 
kinetics 


In  order  to  simplify  the  numerical  calculation,  we  assume 
a  plug  flow  reactor  (PFR)  and  a  one-dimensional  (in  the  ax¬ 
ial  direction)  homogenous  reactor  model  for  the  materials 
and  energy  balances.  The  steady- state  changes  in  the  molar 
flow  rate  of  component  i  ( F* )  and  temperature  (7)  down  the 
length  of  the  PFR  are  calculated  by  numerical  integration  of 
ordinary  differential  equations.  The  independent  variable  is 
the  reactor  catalyst  weight  W,  with  the  limits  of  integration 
between  W=  0  and  W=WCSL t.  The  equations  describing  the 
reformer  can  be  expressed  as  follows: 


W 

dW 


Y  vjiri 

j 


d T  _  Ej(- ^H29,j)rj  +  4 U(T  -  Tw)/pDl 

d W  ~  Ei  FiCw 

where  Ft  is  the  molar  flow  rate  of  component  i  at  the  axial 
position;  W  the  catalyst  weight;  rj  the  rate  of  reaction  7;  pp  the 
stoichiometric  coefficient  of  component  i  for  the  reaction  7;  T 
and  7w  the  reaction  temperature  and  metal  wall  temperature, 
respectively;  A 7/2987  ^eat  °f  reacti°n  for  reaction 7;  U  the 
heat-transfer  coefficient;  D\  the  inner-diameter  of  the  reactor; 
Cpi  the  heat  capacity  of  component  /;  p  is  the  density  of  the 
catalyst.  This  is  a  heated  PFR  with  heat  transferred  from  the 
metal  wall  to  the  reactor.  Because  an  outer  coil  is  placed  to 
heat  up  the  reactant,  the  temperature  of  the  reactor  wall  is 


Table  1 

Rate  expression  for  reactions  in  reformer 

Reactions  Kinetics 


r\  (mol  (kg-cat  min)  1 )  CH4  +  H2O  CO  +  3H2  r\  =  - 

p2.5  / 

h2  \ 

r2  (mol  (kg-cat  min)  - 1 )  CO  +  H2O  CO2  +  H2  r2  =  - 

Ph2  ( 

rs  (mol  (kg-cat  min)  - 1 )  CH4  +  2O2  =>  CO2  +  2H2  r2  =  - 

(1 +K 


modelled  as: 

^condx^  =  4U(T  -  Tw)/pDl  +  4 U(TA  -  Tw)/pDQ  (3) 
d  W 

where  kcond  is  the  heat  conductivity  of  the  reactor  wall;  7a  the 
temperature  of  the  pipe  that  is  coiled  down  the  outside  of  the 
reactor,  and  Do  is  the  outer-diameter  of  the  reformer.  Note 
that,  similar  models  (Eqs.  (1)  and  (2))  are  employed  for  the 
water-gas  shift  reactors  and  preferential  oxidation  reactors, 
but  adiabatic  operation  is  assumed. 

The  reactions  taking  place  in  the  reformer  are: 

ri  :CH4  +  H2O^CO  +  3H2, 

A/7298  =  206  (kJ  mol-1)  (4) 

r2  :  CO  +  H20  <£>  C02  +  H2, 

A//298  =  — 41.2(kJmol—1)  (5) 


r3  :  CH4  +  2H2  <*■  C02  +  2H20, 

A//298  =  -8  1  0  (kJ mol-1)  (6) 


Here,  it  is  assumed  that  carbon  formation  and  carbon  reform¬ 
ing  can  be  neglected  [10].  For  methane  reforming,  the  kinetic 
expression  of  Xu  and  Froment  [15-17]  is  used,  Eqs.  (4)  and 

(5) .  For  methane  combustion  over  a  supported  Pt  catalyst  (Eq. 

(6) ),  the  rate  expression  of  Trimm  and  Lam  [18]  is  adopted 
here.  The  kinetics  of  these  reactions  are  listed  in  Table  1. 

The  reaction  that  takes  place  in  the  water-gas  shift  reactors 
is  the  same  as  Eq.  (5)  except  that  a  different  type  of  catalyst 
is  used.  The  rate  expression  of  Choi  and  Stenger  [9]  is  used 
to  model  the  reactor,  i.e., 


rws  =  ^ws 


f  PcqPy32o  ^h2  ^co2 

V 


Generally,  the  carbon  monoxide  concentration  out  of  the 
water-gas  shift  reactors  is  still  too  high,  i.e.,  at  a  level  that 
can  be  poisonous  to  the  precious  metal  catalyst  in  the  anode 
of  the  PEMFC.  Typically,  preferential  oxidation  is  used  to 
reduce  the  CO  concentration  down  to  the  ppm  level.  In  the 
PROX  reaction,  small  amount  of  air  (oxygen)  is  introduced 
to  oxidize  CO  directly  to  form  carbon  dioxide.  Note  that,  in 
addition  to  the  oxidation  of  CO,  hydrogen  is  also  oxidized  in 
the  same  reactor: 


co  +  Jo2^co2, 


AH298K  =  -283 (kJ mol-1)  (8) 


Sources 


k  1  /JCH4  /JH2 O  -  K\  PC O  Pf 


Ph2o 


1+kcoPco+kh2  pn2  +  ^ch4  ^ch4 +^h2o 


h  Pc O  Ph2 O  -  Pc 02  Ph2  / ^eq 


/  \  2 
1  +  ^CO^CO  +  Vh2  ^h2+^ch4  ^ch4+^h2o 

k]aPcU4  Pq2  hbPCH4Po2 


4 


OX  ,  rO XnA  7 

CH4  Pqh4  +K02  °2  > 


+ 


(l+K°fPCH.+K/5P00 


OX 


ch4 


02 


Xu  and  Formet  [15] 


Xu  and  Formet  [15] 


Trimm  and  Lam  [18] 
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Fig.  5.  Schematic  of  fuel  processor  for  simulation. 


H2  +  -O2  =>  CO2,  A//298K  =  — 243  (kJ  mol  1 )  (9) 

Only  few  papers  address  the  simultaneous  H2  oxidation 
along  with  CO  oxidation  for  the  PROX  reaction.  In  this  work, 
a  constant  selectivity  is  assumed  [11].  The  rate  expression 
of  Amphlett  et  al.  [19]  is  used  for  CO  oxidation  and  40% 
selectivity  is  also  assumed.  Thus: 

r  prox,  1  =  ^prox,  1  PcO  (10) 

^prox,2  =  1.5  X  T prox?  1  (11) 

3.3.  Cooling  device  models 

The  reactor  temperatures  in  the  fuel  processor  are  arranged 
in  the  descending  order,  so  heat-exchanger  and  cooling  de¬ 
vices  are  needed.  Therefore,  a  heat-exchanger  is  installed 
between  the  reformer  and  the  HTS 1  and  the  syngas  out  of  the 
reformer  exchanges  heat  with  the  fresh  liquid  water  feed,  as 
shown  schematically  in  Fig.  5.  This  is  a  standard  exchanger 
and  the  exit  temperature  of  the  syngas  (rexout)  can  be  com¬ 
puted  directly  from  the  energy  balance.  This  is  followed  by 
direct  injection  of  water  to  cool  down  the  syngas  to  the  desired 
temperature.  The  amount  of  liquid  water  required  (Fu2ojnj) 
for  cooling  to  Thtsi  can  be  found  by  solving  the  following 
energy  balance: 

r  Thtsi 

/•  Fi, ATR,out  /  CpidT 

i  4  Texout 

+FH2o.,>y[4i h2o(373  -  298)  +  AHjf 


r  Thtsi 

+  /  CpHjodr]  =  0  (12) 

where  Cp  HoQ  and  Cp,h2o  are  the  liquid  and  vapour  phase  heat 
capacity  of  water  (Jmol-1  K_1),  respectively,  and  AH^0 

represents  the  latent  heat  of  water  (Jmol-1).  Similarly,  the 
water  injected  between  the  HTS  1  and  HTS 2,  HTS 2  and  LTS, 
and  LTS  and  PROX  can  also  be  evaluated  using  Eq.  (12). 

3.4.  Data  regression  for  steady -state  model 
3.4.1.  Reforming  unit 

The  experiment  is  carried  out  in  a  sequential  manner.  It 
follows  the  ATR  pathway  initially  and  then  switches  to  the  SR 
reaction  by  cutting  off  the  oxygen  feed  (Table  2).  From  the 
conversion  data  and  the  feed  flow  rate,  the  material  balances 
can  be  obtained  according  to  the  reactions  (r\  —  r^  in  Eqs. 
(4)  to  (6).  The  inlet  and  outlet  flow  rates  for  both  reaction 
pathways  are  listed  in  Table  2. 

Before  addressing  the  details  of  the  data  fitting,  several 
observations  must  be  made.  First,  it  is  necessary  to  start  from 
kinetics  data  regression  for  the  SR  because  only  two  reac¬ 
tions  are  involved  (n  —  r^  in  Eqs.  (4)  and  (5)).  Second,  the 
composition  is  more  important  than  the  temperature  because 
the  heat  loss  is  not  considered  in  the  reactor  models.  Because 
the  burner  is  not  explicitly  modelled,  the  metal  wall  temper¬ 
ature  (Eq.  (3))  in  the  reformer  inlet  (7w(0))is  used  to  infer 
the  effect  of  the  burner;  a  simple  model  can  be  employed. 
In  other  words,  the  inlet  metal  wall  temperature  is  taken  as 
an  adjustable  parameter.  In  addition  to  7w(0),  the  kinetics 
data  from  the  literature  are  also  modified  to  obtain  better  data 


Table  2 

Inlet  and  outlet  compositions  and  temperatures  data  for  ATR  and  SR  reforming 


Reaction 

pathway 

T  (°C) 

CH4a 

(mol  min- 1 ) 

h2o 

(mol  min-1) 

cob 

(mol  miif  1 ) 

co2 

(mol  min-1) 

°2C 

(mol  min- 1 ) 

h2 

(mol  min-1) 

n2 

(mol  min- 1 ) 

CH4  burner 
(mol  min- 1 ) 

ATR  inlet 

667 

0.48 

0.883 

0 

0 

0.22 

0 

0.83 

0.08 

ATR  outlet 

557 

0.0313 

0.504 

0.204 

0.245 

0.011 

1.172 

0.83 

— 

SR  inlet 

867 

0.48 

1.44 

0 

0 

0 

0 

0 

0.18 

SR  outlet 

517 

0.202 

0.973 

0.087 

0.190 

0 

1.022 

0 

— 

a  97%  Conversions  for  ATR  and  86%  conversion  for  SR  (dry  basis). 
b  Outlet  concentration  of  8.2%  for  ATR  and  5.8%  for  SR  (dry  basis). 
c  02  conversion  of  95%. 


48 


S.-T.  Lin  et  al.  /  Journal  of  Power  Sources  148  (2005)  43-53 


Table  3 


Kinetic  parameters  from  literature  and  from  regression 


Pre-exponential  factor 
(rate  constant) 

Literature 

Regression 

r\  forward 

7.02  X  1016 

6.32  x  1016 

r\  reverse 

5.862  x  103 

1.759  x  103 

r2  forward 

3.258  x  107 

2.77  x  106 

r3  forward  a 

4.866  x  107 

1.56  x  108 

r3  forward  b 

4.092  x  107 

1.31  x  108 

fitting.  Typically,  the  pre-exponential  factors  are  adjusted  to 
match  the  experimental  results  [16]  while  keeping  the  ac¬ 
tivation  energies  and  adsorption  constants  unchanged  [10]. 
Hence,  the  data-fitting  algorithm  of  the  reformer  has  the  fol¬ 
lowing  steps: 

(i)  Starting  with  the  SR  pathway,  7w(0)  is  changed  to  match 
the  outlet  temperature  of  SR  and  the  pre-exponential 
factors  of  r\  and  r2  (Eqs.  (4)  and  (5))  are  adjusted  to 
match  the  outlet  component  flow  rates. 

(ii)  While  keeping  the  kinetics  parameters  of  r\  and  r 2  fixed, 
7w(0)  and  the  pre-exponential  factor  of  r 3  (Eq.  (3))  are 
varied  to  match  the  outlet  temperature  and  component 
flows  for  the  ATR  pathway. 

(iii)  7w(0)  is  constrained  by  establishing  a  linear  equation  to 
relate  7w(0)  with  T(0). 

(iv)  The  kinetics  parameters  of  r\  and  r 3  are  modified  so  that 
the  errors  are  of  the  same  order  of  magnitude  for  the  SR 
and  ATR  pathways. 

The  kinetics  parameters  of  all  three  reactions  from  the 
literature  and  from  the  data  regression  are  shown  in  Table  3. 
The  results  indicate  that  most  of  the  regressed  pre-exponential 
factors  are  quite  close  to  the  literature  data  and  have  the  same 
activation  energies  and  absorption  constants.  The  relationship 


between  7w(0)  and  T( 0)  can  be  expressed  in  °C  as: 

7w(0)  =  1.65  x  r(0)  -  864  (13) 

The  temperature  and  component  flow  rates  along  the  ax¬ 
ial  position  for  the  ATR  and  SR  pathways  are  shown  in 
Figs.  6  and  7,  respectively.  The  outlet  flow  conditions 
predicted  are  roughly  the  same  as  the  experimental  data 
(see  open  circles  in  Figs.  6  and  7),  but  the  outlet  tem¬ 
perature  predicted  by  the  model  is  higher  than  the  ex¬ 
perimental  data,  especially  in  ATR  pathway.  This  may  be 
the  result  of  heat  loss  or  an  inappropriate  location  for  the 
thermocouple. 

The  composition  and  temperature  profiles  in  Figs.  6  and  7 
indicate  that  high  inlet  temperature  leads  to  fast  reaction,  es¬ 
pecially  for  r\ .  This  results  in  a  rapid  decrease  in  the  reactor 
temperature  and  a  fast  increase  in  the  yield  of  hydrogen,  as 
can  be  seen  in  Figs.  6  and  7  for  both  ATR  and  SR  pathways. 
For  the  ATR,  the  oxidation  reaction  (7*3)  leads  to  an  increase 
in  the  reactor  temperature,  so  the  yield  of  hydrogen  and  the 
conversion  of  methane  continue  to  improve  towards  the  mid¬ 
dle  of  the  reactor,  as  shown  in  Fig.  6.  On  the  other  hand,  most 
of  the  methane  is  converted  very  early  in  the  PFR  for  the  SR 
pathway  (Fig.  7).  The  concentration  of  carbon  dioxide  in¬ 
creases  rapidly  initially  for  both  ATR  and  SR  pathways.  For 
the  SR  pathway,  a  lower  temperature  towards  the  middle  of 
the  reactor  favours  V2  (Eq.  (5)),  so  the  concentration  of  carbon 
dioxide  decreases  towards  the  reactor  exit.  For  the  ATR  path¬ 
way,  however,  the  concentration  increases  monotonically  as 
a  result  of  the  temperature  increase  towards  the  end  of  the 
PFR. 

3.4.2.  Gas-cleaning  unit 

After  setting  up  the  steady-state  reactor  models  for  the 
reformer  under  different  operations,  a  similar  approach  can  be 
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Fig.  6.  Temperatures  (metal  wall  and  reactor),  CH4  conversion,  CO  concentration  and  molar  flow  rates  over  axial  position  of  ATR  pathway  (open  circles 
indicate  experimental  results). 
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Fig.  7.  Temperatures  (metal  wall  and  reactor),  CH4  conversion,  CO  concentration  and  molar  flow  rates  over  axial  position  of  SR  pathway  (open  circles  indicate 
experimental  results). 

Table  4 


Experimental  temperature,  composition  and  kinetics  parameters  for  water-gas  shift  reactors  for  ATR  pathway 


Reactor 

7m  (°C) 

CO  (%;  dry  basis) 

H2  (%;  dry  basis) 

Pre-exponential  factor  (literature  [9]) 

Pre- exponential  factor  (regression) 

HTS1 

350 

2.2-2. 5 

44-46 

4.943  X  106 

0.989  x  106 

HTS2 

320 

1.5-1. 7 

45-47 

4.943  x  106 

0.989  x  106 

LTS 

260 

1. 1-1.3 

48-50 

4.943  x  106 

1.285  x  106 

applied  to  the  gas-cleaning  unit.  The  HTS1,  HTS2,  LTS  and 
PROX  reactors  are  arranged  in  series,  so  the  regression  can  be 
carried  out  sequentially.  Also  note  that  all  these  reactors  are 
assumed  to  be  adiabatic.  For  the  water-gas  shift  reaction  (Eq. 
(5))  with  different  catalysts,  Table  4  gives  the  experimental 
data  for  the  ATR  pathway.  As  for  the  reformer  example,  the 


pre-exponential  factor  (Eq.  (7))  is  adjusted  for  data-fitting  for 
the  HTS1  and  HTS2  while  keeping  the  activation  energy  and 
equilibrium  constant.  The  procedure  is  repeated  for  the  LTS. 
The  results  in  Table  4  indicate  that  the  pre-exponential  factors 
from  data  regression  are  of  the  same  order  of  magnitude  as 
those  reported  in  the  literature. 
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Fig.  8.  Temperature,  CO  and  H2  concentration  and  molar  flow  rate  over  reformer,  HTS1,  HTS2  and  LTS  reactors  over  axial  position  for  ATR  pathway. 
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Fig.  9.  Temperature,  CO  and  H2  concentration  and  molar  flow  rate  over  reformer,  HTS1,  HTS2,  and  LTS  reactors  over  axial  position  for  SR  pathway. 


Steady-state  simulation  results  for  the  water-gas  shift  re¬ 
actions  are  given  in  Figs.  8  and  9.  As  a  result  of  exothermic 
reaction  (Eq.  (5)),  the  reactor  temperatures  go  up  along  the 
axial  position  for  all  three  reactors.  A  significant  temperature 
rise  is  observed  for  the  HTS 1 ,  because  of  the  large  conversion 
of  carbon  monoxide,  and  this  is  true  for  both  the  ATR  and  the 
SR  pathways.  By  comparison,  the  ATR  pathway  (Fig.  8)  has  a 
larger  temperature  increase  (Fig.  9)  because  of  its  higher  con¬ 
centration  of  carbon  monoxide  in  the  reactor  inlet.  Finally, 
two  PROX  reactors  are  installed  for  the  further  consumption 
of  carbon  monoxide.  Because  the  concentration  out  of  the 
PROX  is  not  quite  stable,  an  average  value  of  40  ppm  has 
been  taken  for  the  ATR  pathway  to  back-calculate  the  kinet¬ 
ics  parameters.  This  completes  the  steady-state  modelling  of 
the  entire  fuel  processor. 


4.  Sensitivity  analysis  and  optimization 

Once  the  steady-state  model  is  available,  it  is  possible 
to  predict  the  behaviour  of  the  experimental  fuel  processor. 
Design  parameters  such  as  process  configuration,  catalyst  se¬ 
lection,  reactor  design  and  operating  parameters  will  affect 
the  performance  of  a  fuel  processor.  Here,  we  are  interested 
in  the  operational  aspect  of  the  experimental  fuel  processor 
and  the  operating  parameters  studied  include  the  following: 


(i)  The  inlet  temperature  of  the  reformer,  Tm,  is  an  im¬ 
portant  operating  parameter.  A  high  Tm  favours  the  en¬ 
dothermic  steam  reforming  reaction  (Eq.  (4))  and  this 
will  result  in  an  increase  in  the  yield  of  hydrogen.  There 
are  disadvantages,  however,  in  the  form  of  higher  en¬ 
ergy  input  and  higher  CO  concentration  at  the  reactor 
effluent. 


(ii)  The  kinetics  in  Eqs.  (4)  and  (5)  indicate  that  a  high  water- 
to-carbon  feed  ratio  favours  H2  generation  reactions,  but 
it  also  requires  more  heat  input  from  the  burner. 

(iii)  An  increase  in  the  oxygen-to-carbon  feed  ratio  favours 
the  combustion  reaction,  Eq.  (6),  so  it  will  lead  to  a 
higher  reactor  temperature  and  the  endothermic  steam 
reforming  reaction  (Eq.  (4))  is  favoured.  Due  to  a  larger 
share  of  the  combustion  reaction,  less  heat  input  is 
needed. 


Efficiency  is  generally  used  to  evaluate  the  performance 
of  a  fuel  processor  [1,11,12],  and  is  defined  as: 


Efficiency  (EFF%)  = 


LHV  ofH2 

Heat  input  +  LHV  of  CH4 


(14) 


where  LHV  is  the  lower  heating  value  and  ‘heat  input’  is  the 
heat  supplied  from  the  burner  to  the  reformer. 


4.1.  Effect  of  reformer  inlet  temperature  Tin 

The  effect  of  Tm  on  the  experimental  fuel  processor  with 
the  ATR  and  SR  pathways,  while  keeping  other  operating 
parameters  constant,  is  shown  in  Fig.  10.  The  arrows  in  Fig.  10 
indicate  the  base  case  (experimental)  operating  point.  As  Tm 
increases,  the  endothermic  reaction  (Eq.  (4))  is  favoured.  This 
leads  to  better  methane  conversion  and  a  higher  hydrogen 
production  rate  in  both  the  ATR  and  SR  pathways,  as  shown 
in  Fig.  10.  At  the  same  time,  more  CO  is  also  produced  (Eq. 
(4))  and  this  results  in  a  heavier  load  on  the  gas-cleaning 
unit.  The  data  in  Fig.  10  also  reveal  that  the  SR  pathway  is 
more  sensitive  to  changes  in  inlet  temperature  than  the  ATR 
pathway.  As  for  efficiency,  a  higher  Tm  implies  a  larger  heat 
input  along  with  a  higher  rate  of  hydrogen  production.  For 
the  ATR  pathway,  a  rather  flat  efficiency  profile  is  observed 
and  the  optimum  value  is  located  at  7in  =  750°C.  The  SR 
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Fig.  10.  Effects  of  reformer  inlet  temperature  on  CH4  conversion,  H2  production  rate,  CO  concentration  and  efficiency  of  reformer  for:  (A)  ATR  pathway  and 
(B)  SR  pathway. 


pathway,  on  the  other  hand,  shows  a  significant  increase  in 
the  efficiency  as  the  result  of  a  drastic  increase  in  the  hydrogen 
production  rate. 

4.2.  Effects  of  water-to-carbon  feed  ratio 

The  effect  of  the  H20:CH4  ratio  is  explored  with  a  con¬ 
stant  CH4  feed  rate,  as  shown  in  Fig.  11.  For  the  ATR 
pathway,  an  increase  in  the  ratio  favours  the  reforming  and 


water-gas  shift  reactions  (Eqs.  (4)  and  (5))  because  of  ex¬ 
cess  water  (in  a  relative  sense).  This  leads  to  an  initial  in¬ 
crease  (H20:CH4  <2.5)  in  the  hydrogen  yield  followed  by 
a  decrease  (F^ChCFLt  >2.5).  The  reason  is  that  a  further  in¬ 
crease  in  the  F^ChCFLt  ratio  suppresses  the  combustion  re¬ 
action  (Eq.  (6))  and,  subsequently,  lowers  the  reactor  tem¬ 
perature.  Thus,  an  almost  decreasing  trend  in  the  efficiency 
is  also  observed  as  the  ratio  increases  (the  optimum  is  lo¬ 
cated  at  H20:CH4  =  1.8).  The  SR  pathway,  on  the  other  hand, 


Fig.  11.  Effects  of  FbOiCtT  feed  ratio  on  CH4  conversion,  H2  production  rate,  CO  concentration  and  efficiency  of  reformer  for:  (A)  ATR  pathway  and  (B) 
SR  pathway. 
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(a)  (b)  (c) 


Fig.  12.  Effects  of  02:CH4  feed  ratio  on  CH4  conversion,  H2  production  rate,  CO  concentration,  energy  supplied  and  efficiency  of  reformer  for  ATR  pathway. 


shows  a  more  favourable  trend  for  a  high  H20:CH4  ratio.  The 
reason  is  quite  clear,  namely  the  combustion  reaction  is  not 
involved. 

4.3.  Effect  of  oxygen-to-carbon  feed  ratio 

The  C>2:CH4  ratio  is  also  an  important  factor  for  the  ATR 
pathway.  The  influence  of  this  ratio  with  a  constant  CH4 
feed  flow  rate.  As  the  ratio  increases,  the  rate  of  the  com¬ 
bustion  reaction  (Eq.  (6))  increases  and  consequently  heats 
up  the  reactor.  This  leads  to  a  higher  hydrogen  generation 
(Eq.  (4))  and  a  higher  CO  concentration.  As  a  result  of  the 
higher  yield  of  hydrogen,  the  efficiency  goes  up  as  02:CH4 
ratio  increases  and  it  reaches  a  flat  maximum  at  a  ratio  of 
0.45-0.6. 

4.4.  Optimization  of  operating  parameters 

Before  determining  the  optimization  stage,  the  base  case 
operating  condition  (i.e.,  the  arrows  in  Figs.  10-12)  is  exam¬ 
ined.  The  ATR  pathway  has  a  H2  yield  of  1.28  mol  min-1, 
a  CO  concentration  of  40  ppm  and  an  efficiency  of  66%.  By 
contrast,  the  SR  base  case  gives  a  H2  yield  of  1 . 1 8  mol  min- 1 , 
a  CO  concentration  of  30  ppm  and  an  efficiency  of  51%. 
This  corresponds  to  a  power  generation  of  2-2.2  kW.  Be¬ 
cause  the  ATR  reaction  has  a  higher  efficiency  and  a 
larger  hydrogen  yield  with  a  lower  inlet  temperature,  this 
pathway  is  chosen  for  subsequent  optimization.  Also  note 
that,  actually,  the  base  case  operating  condition  of  the 
ATR  is  quite  close  to  the  optimum,  as  can  be  seen  from 
Figs.  10-12. 

Because  the  aim  is  to  find  a  long-term  operating  con¬ 
dition  with  less  expensive  construction  materials,  several 


constraints  are  imposed.  First,  the  reformer  inlet  temper¬ 
ature  is  restricted  to  717  °C  (a  bit  lower  than  720  °C)  to 
prevent  unnecessary  breakdown  over  a  long  period  of  op¬ 
eration  (the  construction  material  for  the  reformer  is  stain¬ 
less  steel  304).  Second,  the  CO  concentration  in  the  syn¬ 
gas  out  of  the  fuel  processor  is  limited  to  45  ppm.  Once 
these  constraints  are  in  place,  the  optimization  variables  in¬ 
clude  02:CH4  feed  ratio,  H20:CH4  feed  ratio  and  Tm.  The 
objective  is  to  maximize  the  efficiency  by  adjusting  these 
three  optimization  variables  simultaneously.  The  results  show 
that  the  optimum  operating  conditions  are:  H20:CH4  =  1.45, 
02:CH4  =  0.45  and  7jn  =  7 17  °C  with  a  hydrogen  production 
rate  of  1.38  mol  min-1.  This  corresponds  to  an  efficiency  of 
68.4%.  This  constrained  optimization  procedure  leads  to  a 
slight  improvement  in  the  efficiency  of  the  fuel  processor, 
viz.,  2.5%.  This  is  because  the  base  case  operating  condition 
of  the  ATR  is  quite  close  to  the  optimum,  as  can  be  seen  in 
Figs.  10-12. 

5.  Conclusions 

Experimental  runs  have  been  performed  on  a  methane  fuel 
processor  from  start-up,  to  autothermal  reforming,  then  to 
steam  reforming.  Next,  steady- state  models  have  been  estab¬ 
lished  to  describe  the  reforming  and  gas-cleaning  units.  A 
systematic  procedure  is  devised  to  adjust  the  kinetic  param¬ 
eters  such  that  a  reasonable  behaviour  description  can  be  ob¬ 
tained.  Based  on  the  model,  sensitivity  analysis  is  performed 
and  the  optimum  operating  condition  can  be  determined.  The 
results  of  the  constrained  optimization  indicate  that  further 
improvement  can  be  achieved  for  the  experimental  fuel  pro¬ 
cessor. 
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